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The present work reports a general method for the calculation of the polarizability of a trun- 
^ ' cated sphere on a substrate. A multipole expansion is used, where the multipoles are not necessarily 

^ , localized in the center of the sphere but can freely move on the revolution axis. From the weak for- 

mulation of the boundary conditions, an infinite set of linear equations for the multipole coefficients 
is derived. To obtain this set, the interaction between the island and the substrate is taken into 
account by the technique of image multipoles. For numerical implementation, this set is truncated 
at an arbitrary multipole order. The accuracy of the method is jugded through the stability of the 
jrt ', truncated sphere polarizability and the the fulfilment of the boundary conditions which are demon- 

strated to be satisfied in large regions of the parameter space. This method brings an improvement 
with respect to the Bedeaux's case MM where the multipoles are located in the center of the sphere. 

a 

O ' PACS numbers: 



> . I. INTRODUCTION 

0^ 



The in situ characterization of the growth mode of a thin film in vacuum conditions is a long pending question. 

\^ , Indeed, a doubt is always shed on the ex situ observations of deposits, often performed by electron microscopy, which 

f^ ' raay not be representative of the studied systems during their growth. Several diffraction techniques are suited to 

0^ ', examine the structure of the growing films. For those many techniques probing surfaces by means of charged species, 

On ' such as low energy electron diffraction and refiection high energy electron diffraction, the charge build-up prevents 

a straightforward interpretation of the results in the case of insulating substrates. Only electromagnetic probes or 

neutral atoms diffraction can analyse most systems in a non disturbing way. However, neutral atoms are only sensitive 

to the top most layer of a substrate while, despite its impressive success in the field of surface science, the grazing 

X-ray diffraction technique [0 is limited in its application because it implies the use of high intensity synchrotron 

2 ' radiations sources. Universally used in the field of thin film growth, near field microscopies as atomic force and 

Q ] scanning tunneling microscopies may imply a perturbation of the growth and are hampered by the tip artefact. 

O ■ Moreover, the rather long time needed to record an image with these methods often prevent them to be used as a real 

^ I time techniques. These limitations likely explain the development of analytical methods based on optical probes for 

monitoring the growth of a film and determining the thickness of the deposited layers. Indeed, UV-visible tools fulfill 

the in situ and non disturbing conditions for the examination of deposits during their formation, combined with a 

simple and versatile use on most substrates. In particular, optical methods, which have also the capability of being 

run in any environment, are widely spread in semiconductor technology Q. 
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At the thermodynamical equilibrium, a film wetting the substrate on which it is deposited is expected to grow 
layer by layer, following a so-called Frank- Van der Merwe growth mode. On the contrary, when a substrate is poorly 
wetted by the deposited material, three-dimensional clusters are formed in a Volmer- Weber growth mode. In an 
intermediate case, the Stransky-Krastanov growth mode, the substrate is covered by a few wetting monolayers before 
the adlayer relaxes to give rise to clusters formation. In addition, in many cases, the morphology of deposited films not 
only depends on thermodynamics, via the surface and interface energies, but also on kinetics through energy barriers, 
diffusion coefficients and intensities of the impiging fluxes. At temperatures and fluxes such that the supersaturation is 
very high, growth of thin films can give rise to the formation of clusters whose shapes are very far from those expected 
in equilibrium conditions. Therefore, it is often very hard to distinguish between the different growth modes at the 
earliest stages of the formation of a film. The knowledge of the cluster mean diameter, aspect ratio, and density during 
deposition would contribute a breakthrough in the understanding of the growth. In this context, surface differential 
reflectivity has been proved to be a powerful tool for the determination of the shape ratio [g|-0| and of the size of the 
metallic clusters [Blp|, even within the monolayer range. However, a quantitative analysis of the optical spectra is still 
at its beginning [g|. 

Since the pioneering work of Maxwell Garnett ||l^] at the turn of the century, there has been large scientific interest 
in the optic al p roperties of metallic clusters |0] . Their optical behaviours are driven to a large extent by the Mie 
excitations |12f| which can be viewed as surface plasmon-polaritons. If, for isolated clusters with simple shapes, 
such as spheres or spheroids in vacuum, the exact solution of the Maxwell equation is known, the difficulty of a 
reliable description of the optical properties of particles dramatically increases for interacting particles of complex 
shapes, either in a matrix or on a surface. Eventhough the Maxwell Garnett effective medium theory and other such 
theories [|3[ have been quite successful in tackling these questions, an accurate description of the macroscopic optical 
properties of the particles, as absorption and reflection, requires a more sophisticated approach than the mean field 
theory to find the renormalized polarizability which governs the far field behaviour. Moreover, for clusters deposited on 
a surface, a quantitative description of the optical properties of the thin film is not only hampered by the interactions 
between aggregates but also by the mutual interactions between the latter and the substrate 0,0 . 

The cluster-substrate and cluster-cluster interactions can be modelled at the dipolar order as in the theory which 
has been developped by Yamaguchi et al. [nqJlTl . The limitations of such an approach, which does not allow for a 
quantitative description of the size and shape ratio of the clusters of a deposit 0] , prompted Bedeaux and Vlieger |lq ] 
to propose a theory, hereafter referred to as the Bedeaux- Vlieger model, to account for the optical behaviour of a 
granular thin film separating two bulk media. This was done by introducing some auxiliary fields named excess fields 
which enabled to describe the macroscopic optical effect of the film, with no need of a detailed description of the 
spatial variations of these quantities as one moves away from boundary layer |18|| . The optical properties of the film 
is essentially related to four surface susceptibilities 113] which are nothing more than the total integrated excess field. 
These surface susceptibilities which are linked to the island polarizabilities for a discontinuous film govern the Fresnel 
optical coefficients pj|-p2| . The modelling of the dipoles arising in clusters under the light excitation relies on the 
assumption that the mean size of clusters is small compared to the wavelength of the light. Under these conditions, 
retardation effects, which are related to the finiteness of the speed of light, can be neglected around and inside the 
islands. Hence, in the electrostatic limit, the Maxwell equations for non-magnetic materials and the Laplace equation 
for the electric potential describe the same physics and are thus equivalent. Similar ideas had been put forward earlier 
by Kretschmann p3[ , but using a somewhat different formalism. 

This method has been applied to truncated spheres on a substrate by Wind et al. 1^,^]. The Laplace equation 
is solved by means of a multipole expansion technique Q which consists in placing multipoles at the centre of the 
sphere. The interaction with the substrate is taken into account by introducing image multipoles located at the image 
point of the multipoles with respect to the surface. The polarizability of the individual islands is then renormalized by 
taking in account the interactions with the neighbouring islands in a self-consistent approximation pl[ | . This relatively 
simple model works quite nicely under many circumstances PJ22J23| and may reproduce the main experimental optical 
features to a large extent (see Sec. || below). 

However, while trying to fit optical data corresponding to disk-shaped clusters with aspect ratio higher than 2, 
which implies a representation of the clusters by truncated spheres with centres below the surface of the substrate m , 
the numerical computation turns out to be unstable. Beyond the necessity of a general model that could be applied 
to any case, there appears a need for clear criteria to establish the consistency of the calculation. In the next section, 
the example of the optical response of a thin silver film deposited on a magnesium oxide substrate allows to describe 
in more details the position of the problem. 



II. DIFFERENTIAL REFLECTIVITY ON AG/MGO(100) FILMS 

In a previous paper, it has been reported that UV-visible differential reflectance data recorded during the deposit of 
a silver film on a MgO(IOO) substrate in vacuum lead to estimates of the size, aspect ratio and density of the clusters 
film which were in excellent agreement with the values derived from an ex situ scanning electron microscopy study of 
the same sample Q. 

A. Experimental setup and results 

A silver film of average thickness 2 nm was deposited on a single-crystal MgO(lOO) held at T = 300°C under ultra 
high vacuum conditions by Knudsen evaporation. During evaporation, p-polarized light (i.e. with the electric field 
in the plane of incidence) in the UV-visible energy range 1-5 eV is impinged at an incident angle of Oq — 45°. the 
reflected light is recorded in the specular direction. The specular reflectivity, R{lu), can be obtained as function of 
the frequency, lu of the incident light. The recorded relevant quantity is the differential reflection coefficient which is 
defined as 

ARju;) ^ Rju;) ~ Rpju;) 
R{uj) Rf{uj) ' ^ ' ' 

where Rpi^^) is the Fresnel reflectivity of the bare substrate. More detailed descriptions of the experimental setup 
and technique used have been already published |||,g| . 

An ex situ image of the silver film has been collected by means of a field electron gun - scanning electron microscope 
(FEG-SEM). It is shown in Fig. ^a). As can be seen from this image, the deposit consists of small islands of a linear 
size of around 10 nm distributed over the surface of the substrate. Such a Volmer- Weber growth mode is expected 
for the Ag/MgO system since noble metals poorly wet wide band gap oxide surfaces [E6[. 

The experimental differential reflectivity spectra at the end of growth is plotted in Fig. Qa as function of the photon 
energy, E = hu, for an incident angle of 45°. The two pronounced resonances seen in this graph are excited by 
the two components of the p-polarized electric field. The low and high-energy peaks are related to longitudinal and 
transversal plasma oscillations inside the clusters, respectively uM- Their positions in energy are mainly governed by 
the aspect-ratio of the islands, which is given by the ratio of height over lateral size, the electromagnetic coupling of 
the film with the substrate and the interactions between particles 



B. Limits of the simulation 

The solid curve in Fig. [ip represents a simulation result obtained by means of the model [^ derived from the method 
of Wind et al mM . Within the framework of this theory, deposited islands are modelled by a set of identical truncated 
spheres which are placed on a regular square array to simplify the calculation since it has earlier been shown that, 
for low cluster coverages, the optical response of a regular array of clusters with, for example, a square or triangular 
lattice, marginally differs from that of randomly distributed aggregates JM • Moreover, for such a low cluster density, 
higher order interactions between islands are expected to be negligible ]27|. The values for the dielectric functions are 
taken from the literature ^^ with, in the case of silver, a correction to account for the finite size of the clusters p|jll|] 
which reduces the electron mean free path. The island polarizability is evaluated in the quasistatic approximation, 
the use of which is justified by the experimental energy range (1-5 eV, i.e. wavelength in the range 200-1200 nm), by 
a multipolar development of the potential that is truncated at an order M for numerical reasons. This polarizability 
is then renormalized by the inter-island coupling which is accounted for only at dipolar order. The model spectrum 
is chosen by means of a try and error method, the main emphasis being put on trying to reproduce the location of 
the low- and high-energy resonances and their intensity. 

The theoretical result shown in Fig. |l|b corresponds to islands represented by spheres of radius i? = 6.8 nm with 
the centre of the sphere located at O.lli? above the substrate. The islands are placed on a square lattice with a 
inter-island distance (lattice constant) of 19.6 nm which corresponds to a cluster density of 2.610^"'^ cm~^. These 
values nicely compare to those obtained by scanning electron microscopy (Fig. nk), a cluster radius of 8.2 ± 1.5 nm 
and a cluster density of 1.9 ± 0.5 10^^ cm~^ g. The silver coverage which can be derived from the optics compares 
within few per cent with that obtained from the SEM image. Such quantitative agreement is impressive, in particular 
in view of the relative simplicity of the model on which the simulation is based. 

The best fit between the optical spectrum and the model has been found for an aspect ratio (diameter/height) 
of 1.8 (Fig. |l|b). Numerical results are shown in Fig. g as a function of the multipole order within this particular 



geometry. A good convergence is obtained upon increasing the multipole order. The same is observed for any case in 
which the aspect ratio of the deposited cluster is lower than 2 (i.e. when the center of the sphere is above the surface of 
the substrate). However, when performing the calculation for truncated spheres with centers below the surface of the 
substrate, the numerical code does not converge any longer. An example is given in Fig. ^ where the model is applied 
to a cluster geometry similar to that used in Fig. except for the aspect ratio which is now set at a value of 2.2. Upon 
including higher and higher multipole orders in the simulations, the position of the low energy resonance is wandering 
and shows no trend of convergence. At this stage, two cases ought to be considered. When the centre of the sphere 
is beyond the mid distance between the surface of the substrate and the bottom of the sphere, the divergence of the 
calculation happens for mathematical reasons since the image point is outside the physical domain or, in other words, 
outside the sphere. In a different way, when the image point is within the sphere, solutions might exist. The lack of 
convergence likely arises from the conditioning of the matrix. Indeed, the problem can not be solved by increasing 
the number of multipoles since, at some stage, the machine accuracy is overflowed. In both cases, a suggestion to find 
out favourable configurations in which the numerical calculation could be performed more successfully, is to move the 



expansion centre. This point will be examined in Sect. Ill 



Nevertheless, neither the convergence of the model itself nor that of the model with the experimental data can 
guarantee the correctness of the solution. These are only indirect proofs. In the model, potentials are expressed 
as a function of expansion coefficients up to the Af th order, although the optical response, which is derived from 
the dipole polarizability of the deposited clusters, only relies on the first order expansion coefficients. Therefore, 
the convergence of the calculated optical response is necessary but not sufficient to validate the expressions of the 
potentials. At variance, the boundary conditions depend on all the expansion coefficients until the highest order, so 
that a definite evidence of the consistency of the model would be gained by their fulfilment. The question will be 



addressed in Sect. IV where the above Ag/MgO case will be used to test the accuracy of the model. 



III. A MODEL FOR EVALUATING ISLAND POLARIZABILITIES 

The cluster geometry used in the present paper, which is derived from that used by Bedeaux et al. is depicted in 
Fig. 0. It consists of a substrate (medium 2) localized in vacuum (medium 1) and covered with clusters (medium 3) 
whose size is assumed to be small with respect to the wavelength of the incident light. The islands are modelled by 
truncated spheres of radius R. The vacuum-substrate interface is located at z — D, with —R < D < R. The plane 
which passes through the centre of the sphere is defined by z = 0. The dimensionless truncation parameter tr 

tr^^, {-l<tr<l), (3.1) 

describes the degree of truncation. For technical reasons which would become apparent below, the part of the spheres 
lying below the interface is introduced as a distinct medium (medium 4) . This latter medium will finally be attributed 
properties identical to the substrate (medium 2). A main spherical coordinate system, S, is defined such that the 
origin of the radial coordinate r coincides with the centre of the sphere (Fig. |j) . The polar and azimuthal angles are 9 
and (j), respectively. The positive z-axis, which is also the axis of revolution of the system, is pointing downwards into 
the substrate. The impiging light of wavelength A scatters at the surface at an incidence angle 9q. To calculate the 
island polarizability, the polarization of the islands by the incident light is modelled by using a multipolar expansion 
and by introducing the image multipoles with respect to the substrate Q]. The centre of the potential expansion is 
chosen along the z-axis at a position fiR in the main coordinate system S, with 

-l<Ai<l. (3.2) 

The image multipole position becomes fiR, where 

p, = 2tr — /i. (3-3) 

It is convenient to introduce two other coordinate systems, S^, and Sp, whose origins are located at the multipole and 
multipole image centres, respectively. The azimutal angle of the plane of incidence of the light is defined as 4>o- It is 
taken equal to zero in all the numerical calculations which are performed herein. 

A. Multipole expansions of the potential 

If the linear dimension of the island is small compared to the wavelength of the incident light, retardation effects in 
the island and around can safely be neglected. If the island material is assumed to be non-magnetic, then the incident 



magnetic field will be unaffected. Hence the main field is the electric field E(r), which is related to the potential, ip{r), 
in the usual way E(r) = — V^(r). When retardation effects can be neglected or, in other words, when the size of the 
island is supposed to be much smaller than the wavelength, this potential must only satisfy the Laplace equation 



V^V-Cr) ^ 0. 



(3.4) 



The appropriate boundary conditions for this potential come from the continuity of the potential itself and of the 
normal components of the displacement field 



V'/(rs) ^tpiiivs), 
ei{uj) 9„V/(rs) = £ii{^) a„V//(rs). 



(3.5a) 
(3.5b) 



Here, ei(ti;) is the frequency dependent dielectric function of medium i, r^ is assumed to be any point on the interface 
and dn denotes its normal derivative. These conditions should be fulfilled for any two media having a common 
interface. 



To solve Eq. (3.4), with the boundary conditions Eqs. ( |3.5| ), it is convenient to use pq] a multipole expansion for 
the potentials. In medimn 1, the potential then takes the form 



1^0 i^a 

Im hn 



(3.6) 



The spherical coordinates (r^,0p,(/)p) and {rp,,Ofi, cjjfi) are referring to the coordinate systems 5^ and 5p, which are 
centered at the points /i and jl, respectively (Fig. m. The m-summation will be understood to range from m = —I to 
I. The spherical harmonics Yi„i{9,(f>) are normalized according to the convention |M 



Y, 



Im \ 



121 + 1 {l~m)l 
Att {l + my. 



Pr (cose) {-!)"' e'""^, 



where P™ (cos 9) are the associated Legendre polynomials defined as 



2'/! 



dxl+m 



X = cos f 



for TTi > 0, while for m < 



(i — my. 



(3.7) 



(3.8) 



(3.9) 



In Eq. ( p.6[) , "00(1") stands for the potential corresponding to the incident field Eq which, provided clusters are small 
compared to the wavelength of the light, can be approximated by a homogeneous field 



Eq = Ef) (sin do cos 00 , sin ^o sin 0o i cos 9q ) . 
The potential 00(1") related to this incident field, then takes the following form 

0^0 (r) = —rEo (cos 9 cos ^o + sin 9 cos sin 9q cos 0o + sin 9 sin sin 9o sin 0o ) 



(3.10) 



= -rEn\ —- 



V2 cos 9aYi 



10 l 



% {e-'*«Yn{9,^) ~ e'^°Y,MO,cl^)} 



(3.11) 



with Eq = |Eo|. The spherical coordinates {r,9,(j)) are defined with respect to the coordinate system S whose origin 
is at the centre of the sphere. 

For the potential inside the substrate, the solution of the Laplace equation is chosen as 



i'2 (r) = flo - rEo 



27r r 



V2ai cos^o^K 



3in0o{a2e-*^°rii(0,0)~a3e^'^''yi_i(0,0)} 



i#o 
T.<r.r-'-'Y,m{9f.,, 

Im 



(3.12) 



where ai are constants to be determined . No tice that the terms inside the square brackets are a hnear combination 



of the terms involved in the expression ( 3.10 ) of the potential V'o(r). For the potential inside the cluster, medium 3 



and 4, the following expansions are defined 



1^0 1^0 

Im Im 

V'4(r) = h'^ + Y. ^^rn rlYi^idf., 0^) (3.13b) 

Im 

where bo and fog are constants to be determined. 

B. Determination of the expansion coefficients via boundary conditions 

The various unknown expansion coefficients, a^. Aim, Bim, etc. . . are now calculated in terms of the known param- 
eters of the model by imposing the boundary conditions on the considered geometry. Since a formulation with image 
multipoles is used, the boundary conditions at z ^ D are easily satisfied. Their fulfilment leads to the following 
relations 

ao = Eod \ 1 I cos 9o , 

V£2 J 

£l 

01 = — , 

£2 

02 = as = 1, 

£l +£2 

AT _ ^£]_4, 
£l +£2 

The continuity at the same boundary {z = D), but now inside the sphere, gives 

£3 +£4 

bL = -^^Bim. (3.15c) 

£3 + £4 



(3.14b) 


(3.14c) 


(3.14d) 


(3.14e) 


(3.15a) 


(3.15b) 



An explicit expression for the constant feo is given in Eq. (A12) of Appendix 1^. These results are obtained by taking 

in th( 
mage 

IF 



advantage of both the orthogonality of the spherical harmonics and the fact that on the surface of the substrate 
[z = D), as a consequence of the symmetry of the location of the multipoles and i mage multipoles, the following 



relations hold r^ = r^i, 9^ = n — 9p, and 0^ — (f>p,. Note that the relations ( 3.14 )- (|3.15| ) are independant of the 
location of the multipoles, so that they are similar to those obtained by Wind et al. pj under the assumption that 
the multipoles arc located at the centre of the sphere. Finally, there are only two independant classes of expansion 
coefficients, namely Aim and Bim- 

These coefficients can be derived from the relations that express the fulfilment the boundary conditions at the 
surface of the sphere, where r = R. By multiplying all the terms of those relations by the complex conjugate spherical 
harmonic, [!";,„ (0, (/>)]*, and by integrating the resulting expressions over the surface of the sphere, where one again 
takes advantage the orthogonality of the spherical harmonics, the following (infinite) set of equations is obtained 

00 
J2 [ci;^R-''-^Ai,m + D^,R''-'Bi,m] = Hr, (3.16a) 



00 



n. Aiirn + ^ii'-ri ^I'm 



l' = l 



= jr, (3.16b) 



where I = 1,2,..., and tti = 0, 1. The various matrix elements of the above linear system can all be found in 
appendix ^, where their detailed derivation also is presented. The above linear system has only a non-trivial solution 



when m = and m = ±1 (cf. Eqs. (A4). These values are associated with the components of the uniform incident 
field Eq which are orthogonal and parallel to the substrate, respectively. An additional simplification can be made by 
showing that the equations corresponding to to = —1 and rn — 1, r espec tively, are equivalent. 



To allow for a numerical solution of the infinite linear system ( 3.16 ), an upper cut-off M in I and I' has to be 
introduced. The integrals that define the system of equations (A5) are computed by Gauss-Legendre numerical 
integration. Solutions of the linear system lead to values of the expansion coefficients Aim and Bim, which in turn 
can be directly related to the (dipole) polarizability of the island. Those can be shown to be given by [gj 

a± = —fw^ -r^ (3.17a) 

^■|_Eocos6'o 

'""^^" (3.17b) 



-^ii'o sin 00 exp(— i0o) 



Here, a_L and a\\ are the dipole polarizabilities perpendicular and parallel to the interface of the substrate, respectively. 
Due to the presence of the substrate, these are in general quite different. In the present work, the parameter of interest 
is the differential optical reflectivity /S.R{uj) / R{w) defined in Sec. O. This quantity is evaluated by modified Fresnel 
formulae for reflection where the above given polarizabilities appear through the surface susceptibilities iQ . 

The present section was aimed at introducing a general method for the calculation of the polarizability of truncated 
spheres supported by a substrate. The model of Wind et al. [y|2| is just a special case of this approach with /^ = (see 
Appendix ^ for details), at least for tr>Q. Indeed, to obtain the solution for the case in which the centre of the sphere 
is located below the substrate (i^ < 0), these authors |jl|,|| apply to the {t^ > 0) case a coordinate transformation 
method with a permutation of the dielectric constants. This implies that different potential expansions are used in 
the two cases, while, in the present work, the expansions in use do not depend on the location of the multipoles with 
respect to the surface. Criteria are now needed to know how reliable are these models, special attention being paid 
to cases corresponding to negative values of the truncation ratio tr- 

IV. FULFILMENT OF THE BOUNDARY CONDITIONS 

In Sec. y and in the already published work ||^, the method has been only judged on the basis of the convergence 
of the differential reflectivity curves AR{uj) / R{uj) . Indeed, in the model used herein, the reflectivity only relies on 



the lowest order expansion coefficient Aiq and An (cf. Eqs. (3.17)). Therefore, the occurrence of a convergence of 
the refiectivity curves upon, say, increasing the multipole order M, does not provide any proof of the reliability of 
the potentials since those depend on all the expansion coefficients up to the cut-off order. The potentials given by 



the expansions (3.6), ( 3.12| ) and ( p.l3| ), can in principle be calculated in any point of space provided the multipole 



coefficients, Aim and Bim, are known. Since the multipole expansions are indeed solutions of Eqs. ( |3.4[ ) and ( |3.5| ) at 
any M , it is only the fulfilment of the boundaries conditions at the interface between the various media that has to 
be checked. It has to be stressed that this fulfilment is a more severe test for the calculation than just considering the 
convergence of the differential reflectivity curves. It is not only the self-consistency of the numerical implementation 
but that of the method itself that is controlled. Such a study thus provides a powerful and rigorous tool for justifying 
the quality of previous and present numerical simulations based on a multipole expansion. 

The calculation including the multipoles and their images, the boundary conditions at the cluster-substrate interface 
are obeyed by construction. Thus, the only boundary conditions to be checked are at the surface of the sphere. To 
measure the error in these, two error functions are defined 

V>+(r,)-V;_(r,.) 
maxr, ipQirs) 

ta„^(rs) = (-,11 w • (4.1b) 

maxr, (ec'„Vo(rs)) 

The subscripts -|-(— ) correspond to quantity just outside (inside) the surface of the sphere. The maximum value of 
the incidence potential, '00, at the surface of the sphere, and the corresponding quantity for the normal derivative are 
used as normalisation factors. The impact of the parameters of the model, location of the multipole centre, photon 
energy and multipole order on the accuracy of the numerical calculations is examined in the following by considering 
the values taken by the error functions in various representations of the above mentioned Ag/MgO deposit. 



A. Location of the expansion centre 

In the model put forward in the present paper, the multipole expansion point can move along the symmetry z-axis, 
instead of being at the centre of the sphere as in the Bedeaux-Vlieger model. In Fig. || an attempt is made to 
determine the optimal position /i, for a given truncation ratio and value of M = 20 at an energy oi E = 4. 5 eV . 
The mean spatial absolute error in the boundary conditions (|£^(rs)|) a.nd{\£g^^^{rs)\), which are defined in Eq. (ff.]]), 
are shown in Fig. |5| as function of the truncation ratio. It can be observed that the value of /i which leads to the 
minimum error depends on the truncation ratio. For tr > 0, the Bedeaux-Vlieger limit, where n = 0, provides the 
minimum error whereas, for tr < 0, an improvement is obtained by moving the expansion point inside the physical 
domain or, in other words, by using negative fi values. However, for negative truncation ratios, the matrix system 
become dramatically more and more ill-conditioned, a fact that may be due to intervention of an increasing power of 
distance between the surface of the sphere and the expansion points /i and p,. 

In principle, the infinite expansion of potentials on which the method is based should lead to a unique numerical 
value, wherever the location of the expansion centre is. In the case of numerical calculations, where the expansion is 
truncated at an order M, it is not too surprising to find, for a given value of M, a (weak) dependency of the error 
functions on the position of the expansion point, as can be already stressed in Fig. o. A search for an optimal position 
of the expansion centre can be performed for each value of the truncation ratio (Fig. ||). For tr > (see also Fig. ||), 
it is the centre of the sphere (/x = 0). For tr < 0, it appears that the best choice is to place the /i point close to the 
substrate so that fi = tr, that is to say to superimpose the centre of expansion and its image. Roughly speaking, 
for tr < 0, between one and two order of magnitude in the error function can be earned by the present method with 
respect to the Bedeaux-Vlieger model in which /i = 0. 

B. Boundary conditions at the resonance energy 

When the energy of the incident light is close to a resonance of the system, i.e. in our case, close to either the low 
energy resonance around 2.5 eV or the high energy resonance at 3.7 eV, the boundary conditions are even harder to 
satisfy. This is seen from Fig. where the (spatial) mean of the absolute error in the boundary conditions, {\£^{rs)\) 
and (|fa^^(rs)|), are shown as function of energy E = hoj for truncation ratio tr = 0.1 and tr = —0.1. The dips which 
can be observed around the low energy resonance at E = 2.5 eV are believed to come from numerical artefacts. These 
figures highlight the sizable effect of the dielectric function on the quality of the simulation, since the energy range 
E = 2 eV to E — 3.7 eV, where the resonances are peaking, concentrates the highest values of the error functions. 
This is, however, probably not so surprising in view of the fact that the system is close to singular at a resonance 
point. As the truncation ratio is reduced, the overall error seems to be increased, but it is still the regions around the 
resonances that are associated with the largest errors. Roughly speaking, the errors in the boundary conditions are 
raised by an energy independant factor as the truncation ratio is reduced (Fig. |^). 

C. Multipole order 

By increasing the number of multipoles included in the calculation. A/, the fulfilment of the boundary conditions 
should be improved since the size of the function basis increases. Fig. ^ shows the errors {\£^{rs)\) and (|i?a„i/;(rs)|) 
as function of M for the energy E = 4.5 eV and /i = for two truncation ratios tr — 0.1 and tr = —0.1. These figures 
are representative of the cases where fi — 0. For tr > 0, the error decreases upon increasing M, which corresponds to 
the expected behaviour. However, for tr < 0, the error is not seen to decrease upon increasing M. In this case, the 
error functions are quite sensitive to the value of /i. This is illustrated in Fig. ^ where the error functions are given 
for /x- values of 0. and —0.1. It is observed that moving the expansion point away from the center of the sphere may 
dramatically improve the overall error made in the simulation. Finally, it can be noted that, in all cases, a higher 
limit in M is imposed by the vanishing of the matrix conditionning. 

D. Spatial variations of the error functions 

The evolutions of the potential -0 and of its normal derivative ednijj along the surface of the sphere, in the plane 
of incidence of the light (Fig. |l(]| a and F ig. jlQJ b) and in a plane perpendicular to it (Fig. |lO| c and Fig. [lC|d), are given 
as function of the polar angle 6 in Fig. HO^ for the case corresponding to tr = 0.1, E = 4.5 eV, /i = and M — 16. 
Also shown are the error functions (|f^(rs)|) (Fig. lOi and Fig. ^0|c) and (|£a„^(rs)|) (Fig. |l0| b and Fig. p^) in the 



bottom panels of the graphs. For the same set of parameters, we have also found (results not shown) that the errors 
are slowly varying function of the azimuthal angle 4>. The regions which give rise to the largest errors in the boundary 
conditions are concentrated around the top of the sphere and, to some extent, at the interface between the spherical 
cap and the substrate. The same trends can be observed in the global view of the boundary conditions which is 
presented in Fig. n^ for the same case. This view consists in a projection in the plane of incidence of (|5^, (r^)!) (upper 
panel) and {\£driii{^s)\) The direction of the incidence of the light is indicated by an arrow. The error is roughly of 
the same order of magnitude all over the sphere. 

V. CONCLUSIONS 

A generalised method for the calculation of the polarizability of a truncated sphere by means of a multipolar 
description has been presented. The proposed approach, which is an extension of earlier published works due to 
Bedeaux, Vlieger and co-workers [|l|,H, and which contains this model as a special case, is based on the possibility for 
the multipole expansion centre to freely move along the z-axis of the main coordinate system instead of being located 
at the centre of the sphere. This method allows to perform the calculations for cases in which the numerical code is 
badly conditioned when the centre of the truncated sphere is below the plane of truncation. Let us stress that this 
geometry has an strong importance since it corresponds to clusters whose aspect ratio is higher than 2. 

Neither the convergence of the model itself nor the agreement between the data and the model can guarantee that 
the solution is correct, since the reflectivity curves only depend on the lowest order expansion coefficients. A way to 
assess the quality of the numerical calculation up to the upper order has been introduced which consists in judging 
the accuracy of the models on the basis of the fulfilment of the appropriate boundary conditions. It is shown that, 
provided the cluster aspect ratio is not too high, the multipole expansion method represents a rather accurate tool 
for the determination of mean values of the parameters characterizing the deposited clusters. In all cases however, 
the values found for the error functions allow to estimate the validity of the model. The simulation of the optical 
reflectivity in the UV-visible range allows to determine cluster size, shape ratio and density in a rather accurate 
manner, as illustrated herein by the case of a Ag/MgO(100) deposit. The method offers a promising tool for the in 
situ examination of cluster growth on substrates. 
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APPENDIX A: THE MATRIX ELEMENTS 

In this appendix we give the various matrix elements, and the form of the right hand side, appearing in the matrix 
system (Eqs. ^.16 ). With the particular expansions chosen for the potentials ipi{r) {i = 1,2,3,4) (cf. Sec. Ill), and 



the relations Eqs. (3.14) and ( |3.15 ), the boundary conditions at the surface of the substrate, by construction, are 



automatically fulfilled. However, we still need to satisfy the remaining boundary conditions Eqs.(^^) on the surface 
of the sphere. By taking these boundary conditions, multiplying by the complex conjugate spherical harmonic, 
[YimiO, </')]* I cind integrating over all direc tions , i.e. using the weak formulation of the boundary conditions, one is lead 



to the following matrix system (cf. Eqs. (|3.16| )) 



oo 



R Al'm + DJ^,R Bi/,n 



l' = l 



= HP, (Ala) 

J2 [Fri^R-''-^Ai,,n + Gi;5i?''-^i?i'™] = Jr- (Alb) 



i'=i 



This linear system can be used to determine the expansion coefficients Aim and -B/„i for all allowed values of / and m; 
I = 0, ±1, ±2, . . . and ni — 0, ±1, ±2, . . . , ±m. However, for orthogonality reasons on integration over </>, this system 



is reduced to m = —1,0, 1, since only the (un iform) incident field contains these "quantum numbers" (cf. Eq. (B.ll) 
The matrix-elements of the system (Eqs. Al) are thus given by the following expressions 



/^ni /-ni 



2£l 



K^'Wtr) + ^—^ (-1)''+™ Kimitr) + —:^ {KJ^MiU = 1) - Kl^^WU)} 



£l +^2 



El +£2 



T)ra . /-m 



MJ]JMU) + '-^-^ (-1)''+" MJj^MiU) + -^ {MIPMU = 1) - Ml^MU)} 



£3 +£4 



£3 + £4 



ipm /-ni 



2£l£2 
£1 +£2 



7- 771 



Wu = 1) + ei?^ |iI?[M](i.) + (-i)''+™Li7,[^](t.)} 

£1 + £2 •- J 



and 



Gm /-Ti 



'^'''^ NJ^[lA{tr = 1) +£3^M^ [nJ^WU) + (-l/'+"7V;7?[M](ir.)} 



.£3 + £4 £3 + £4 

Here the following notation has been introduced 



/-m 

siv — 



1 [21 + I) {21' + I) [I - m)\ {V - m)\ 



[l + m)\ {V + m)\ 



(A2a) 
(A2b) 
(A2c) 

(A2d) 

(A3) 



and the integrals K, L, M, and N will be defined below. Furthermore, the right-hand side of the system Eqs. (Al) 
is defined as 



Hr = V^ 



bo „ £1 - £2 a 

— - EqU cos6'o 

-K £2 



^ —Eocos9o6o^ 



SoiSom 

£1 -£2 



— hi 

£2 £2 



ivStrCwQioitr) QlQll[tr)f 



EoCOseoSn\e-''^"Si^-e''^°S.ir 



(A4a) 



j A^ JO-TT 

Jr = Y Y^o£i cos 00 SomSu - y Yi?o£2 sin ^q -51, [e-^*° 5i,n - e^*° ^_i^] 



and 



— [(£1 - £2)£;osin0o] [e-"^°CuQn{tr)Sim - e'^\i^'Qi,\tr)5. 
The above equations depend on several types of integrals. They are defined as 



K]?[v]{u) = I dxprix)prl ^ y^-M [xMx,i)]- 



LTl'Mtr) 



dx 



-1 



Pr (X) dr \ Pr I 5I_^ ) [^[^] (,, ,)] ^ 



MiifiviWr) = rdxpr{x)pr ( 4=4= ) [xW(a;,ir , 

-1 VvxW(a;,l)/ 



and 



A^i™M(tr)= / ^dx 



pr{x)dApii 



-ri{tr 



\/x['n\{x,r) 



[xW\{x,r)]' 



(A4b) 

(A5a) 
(A5b) 
(A5c) 
(A5d) 

(A5e) 



All the above integrals, except the Q-integral, have a functional dependency on the function rjitr), which is used as 
a generic case for either /i(tr) or fl{tr), i.e. for the z-coordinate of the position of the multipoles or image multipoles 
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in the main coordi nate s ystem S with the origin in the center of the sphere. Moreover, the functional, x['?](^j'') 
appearing in Eqs. ( A5b )-( A5e| ), is defined as 



X[v]{x,r) =r'^ -2'q{tr)xr + ff{tr), 



cost 



(A6) 



This functional gives the distance to a given point from the location of the multipoles (77 = fi) or image multipoles (77 = 
p) in terms of the coordinates, (r, 9, (f), of the main coordinate system S. Notice that xMCa^j i") is independent of the 
azimuthal angle 0. This is, of course, a consequence of the rotational symmetry around the z-axis of the considered 
geometry. 



We would like to stress that with Eq. (3.9) and the definition of C;™, Eq. (A3), the matrix elements constituting the 
left-hand-side of our linear system Eqs. (A2), do not depend on the actual sign of to, but only its value, i.e. C;™ = C'^y,™ 
with similar expression for the other (left hand-side) matrix elements. For the right-hand-side matrix elements, it 
follows directly from the definitions of these quantities, Eqs. (A4), that they vanish identically for to, 7^ 0, ±1 and 



that they furthermore satisfy the relations Hi exp(i(/ 



--ff; exp(i0o) and J} exp( 



-J; ^exp(i0o)- This 



is so because the right-hand-side of the linear system stands entirely for the incident field, which is the same in the 
present case and the one considered in Ref. |y], and therefore is independant of the location of the multipoles. As a 
consequence it follows that the multipole coefficients fuUfil the following relations 






-icjio 



(A7) 
(A8) 



Hence, the independent multipole coefficients are chosen to those corresponding to to = 0,1. All other coefficients 
either vanish or can be expressed in terms of these. Furthermore, it can be shown that if the physical system does not 
contain any free charges, then the equations with Z = is useless. Thus the linear system Eqs. (Al) will only be of 
interest when Z = 1, 2, 3, . . . and to, = 0, 1. The above findings look like those found in the Bedeaux-Vlieger model [Q. 
We may take the "Bedeaux-Vlieger" limit of our results by placing the multipoles in the center of the sphere, i.e. 

(A9) 
(AlO) 



(Alia) 
(AUb) 
(Allc) 
(Alld) 



Using the notation of Wind et al. [pip 






M(ir) 


= 0, 






/i(t.) 


= 2tr. 






e introduced integrals then take the folk 






MJ^WU)-- 


-Ql!.(tr), 






Mipmu)-- 


-TlPiU), 






NIIiMtr) = 


'-I'Qu'iU), 






NIPMitr) - 


- dt TiY^{tr)\t=^ 



The QJ^-integral, which is independent of the position of the multipoles, is the same in both works. The linear system, 
of Ref. ||l|, used to determine the multipole expa nsion coefficients, can now be obtained from our formulae Eqs. ( Al) 
by taking advantage of the relations Eqs. f|A ll|). In this limit our formulae for the linear system reduce, as they 
should, to those presented by Wind et al. IM 



In order to evaluate the potentials in the check of the boundary conditions, the constant bo, appearing in Eq. (B.13a), 
must be determined. This constant is obtained in the same way as the above matrix elements, but here only terms 
containing Yqq{9,(/)) will contribute. The result is 



bo 



£1 
£2 



:Eo{--1] cos0oCoiQoi(ir-) + REoU ( ^ - 1 ) COS0O [1 - CooQoo(^r)] 



£1 

£2 



^ 00 



^ oc 

— = Y^ BuoR^ CI 



K^oi'Mtr) + {-iy":^-^K'[mr) 



£1 +£2 



MS,[^,]iu) + (-iy"^^M^Amr 



2£i 



£3 + £4 



£1 + £2 
2£3 



£3 + £4 



{K",MW-KUl^]itr)} 



{M°,[M](l)-il4,[M](i.)} 



(A12) 
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FIG. 1. Experimental data for a 2 nm thick silver deposit on a magnesium oxide MgO(lOO) substrate: (a) A FEG-SEM image 
of the film. The horizontal bar indicated in the figure corresponds to 200 nm; (b) The experimental differential reflectivity 
spectrum obtained at an incident angle of 45° (circular symbols). A simulation result, using the method of Ref. pJJ3|, is indicated 
by a solid line, and demonstrate good agreement with the experimental data. In the simulation, the truncated spherical particles, 
of radius R = 6.8 nm, were placed on a regular grid of lattice constant 19.6 nm, and the truncation parameter (see text) was 
tr = 0.11 (corresponding to a aspect ratio of 1.80). Notice that the main spectral features and the numerical values of both 
the low and high-energy resonances are well predicted by the theoretical model. 

FIG. 2. Simulated differential reflectivity curves AR{u)) / R{u)) as a function of energy E = hu for differents number of 
multipoles included in the simulation. Clusters are defined by the same numerical parameters as in Fig. hi. 

FIG. 3. Simulated differential reflectivity by the Bedeaux-Vlieger truncated spherical model for the same islands parameters 
as in Fig. hi but for a negative truncation ratio tr — —0.11 (aspect ratio of 2.25). The various curves correspond to different 
choices for the multipole order as indicated in the legend of the figure. Notice that in the region around the low energy resonance 
no convergence seems to be reached by increasing the number of multipoles included in the calculation. The inset shows the 
details of the differential reflectivity curves around the low energy resonance. 

FIG. 4. The cross section of the geometry considered in the present work. Here nR indicates the z-coordinate of the location 
of the multipoles, while fiR is the same quantity, but for the image multipoles. Notice that the multipoles and thus also the 
image multipoles, are always located on the z-axis which is the axis of revolution. The position of the substrate is parallel to 
the xy-plane, and it is located nt z — D, where D is a (signed) real constant. The dielectric functions of the various regions 
are those indicated in the figure by ei{u>), with i = 1, . . . , 4. 

FIG. 5. The mean errors as defind in the text, (|i£^(rs)|) and {|£g,,^(rs)|), in the boundary conditions of (a) the potential 
and (b) its normal derivative, as function of tr for different /i values. The model clusters are defined with the same parameters 
as in the Ag/MgO(100) case (Fig. |l|), with E = 4.5 eV and M = 20. 

FIG. 6. Fulfilment of the boundary conditions as function of fi for an energy E — 4.5 eV, a multipole order M = 20 and two 
given values of tr, 0.1 and —0.1. Error functions on (a) the potential and (b) the normal derivative. 

FIG. 7. The mean error in the boundary conditions as a function of energy for two values of tr — 0.1 and tr = —0.1, in the 
case defined by /x = and M = 20, for (a) the potential and (b) the normal derivative. 

FIG. 8. The evolution of the error in the boundary conditions as a function of the order in the multipole development in the 
case of which E = 4.5 eV, i-i ~ 0, for two given values of tr, 0.1 and —0.1 ((a) potential (b) the normal derivative). 

FIG. 9. The improvment brought in the convergence of the boundary conditions by moving the expansion point in the case 
tr = —0.1, E = 4.5 eV as a function of the multipole order AI. 

FIG. 10. Evolution of the potential (a and b) and its normal derivative (c and d) along the surface of the sphere. The 
errors in the boundary conditions of the potential on the surface of the sphere, £^{rs) (a and c) and of the normal derivative 
£g^^{rs) (b and d), are shown as function of the spherical coordinate 6 for given value of the azimutal angle cfi. The values used 
for the angle (p correspond to the incident plane {(j> = 0) (a and b) and a plane perpendicular to the incident plane {(j> — it/2) (c 
and d). The energy of the incident light is _B = 4.5 eV, the truncation ratio is tr = 0.1 (and fi = 0.) and the number of 
multipoles is M = 16. 

FIG. 11. Global view of the boundary condition errors (a) (|f^(rs)|) and (b) {\£a„ii{i^s)\) when projected into the plane of 
incidence. The parameters used are E = 4.5 eV, f,- — 0.1, and M = 20. 
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Figure 2 
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Figure 3 
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I. Simonscn, R. Lazzari, J. Jupille and S. Roux 
Optical response of supported particles 



22 



w 




e 







2.0 




1.0 




0.0 




-1.0 




2.0 


1 

o 


0.0 


5- 

ro 
til 


-2.0 




Figure 10 
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